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Keeping in mind the current picture of an accelerating and flat Universe, some specific dynamical models of the 
cosmological term A have been selected for investigating the nature of dark energy. Connecting the free parameters 
of the models with the cosmic matter and vacuum energy density parameters, it is shown that the models are 
equivalent. Using the selected models, the present values of some of the physical parameters have been estimated, 
and a glimpse at the past decelerating universe has also been presented. It is observed that most of these cosmological 
parameters nicely agree with the values suggested by the Type la Supernovae and other experimental data. 
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The observations on supernova by the High- z Super- 
"ch nova Search Team (HZT) and the Supernova Cosmol- 
Q_i'Ogy Project (SCP) [1,2] have revealed that, instead of 
slowing down, the expanding Universe is speeding up. 
An intense search is going on, in both theory and ob- 
servations, to unveil the true nature of this acceleration. 
It is commonly believed by the cosmological community 
that a kind of repulsive force which acts as anti-gravity 
! '"I is responsible for gearing up the Universe some 7 billion 
rS years ago. This hitherto unknown exotic physical entity 
is termed as dark energy. 

Now, there can be many variants of dark energy 
which can be responsible for this accelerated universe, 
and variation in the forms of dark energy also exhibit 
variation in expansion rates in different eras. So, there 
may be more than one candidate which can be stamped 
as dark energy. For example, one may select the so- 
called cosmological constant, introduced and later aban- 
doned by Einstein, as dark energy. But selection of 
the cosmological constant as dark energy faces a seri- 
ous fine-tuning problem which demands that the value 
of A must be 123 orders of magnitude and 55 orders of 
magnitude larger on the Planck scale (T ~ 10 19 GeV) 
and the electroweak scale (T ~ 10 2 GeV), respectively, 
than its presently observed value. Moreover, the matter 
and radiation energy densities of the expanding Universe 
fall off as a~ 3 and a~ 4 , respectively, where a is the 
scale factor of the universe, while A remains constant. 
This poses another disturbing fine-tuning problem. For 
these two reasons, at present A with a dynamical char- 
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acter is preferred over a constant A, especially a time- 
dependent A which has decreased slowly from its large 
initial value to reach its small value at present [3]. A 
scalar field <j) with a potential V ((/>), which is known 
as quintessence and decreases slowly with time, may be 
another candidate for dark energy. Quintessence exerts 
negative pressure and is dynamic in nature (recall that 
Aefiective = SirGp^). However, in the present article we 
have considered some phenomenological models of kine- 
matical A which is assumed to be one of the dark energy 
candidates to account for the accelerating expansion of 
the Universe. 

Among the dynamical models of A which are fre- 
quently used in the literature, we have particularly pre- 
sented here three types, viz., A ~ (a/a) 2 , A ~ a/a and 
A ~ p, where a is the scale factor of the Robertson- 
Walker metric and p is the matter energy density. The 
first type of A -model was proposed from dimensional ar- 
guments by Carvalho et al. [4] and Waga [5] and using 
another type of argument by Lima and Carvalho [6] , and 
it was subsequently taken up by several workers [7-11]. 
Using dimensional arguments, Vishwakarma [12] sug- 
gested the A ~ p model, whereas the second model 
mentioned above was dealt by Arbab [13-15] and Over- 
duin and Cooperstock [3]. 

Now, a key to catch up the nature of dark energy lies 
in w , the equation of state parameter which is nothing 
but the ratio of fluid pressure and energy density of dark 
energy, viz., w = p/p. This parameter w has different 
forms in different models. In the present study, using 
the above three forms of A, general solutions of the 
field equations are obtained under the assumption that 
the Universe is flat. Also, particular solutions, wherever 
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needed, are discussed for the specific cases of matter- 
and radiation-dominated universes related to three spe- 
cific dynamic cosmological terms. It is possible to show 
the equivalence of the three models in terms of the so- 
lutions, obtained by connecting the free parameters a, 
[3 and 7 of these models with Q m and , the mat- 
ter and vacuum energy density parameters of the Uni- 
verse, respectively. This will enable us to establish a 
relationship between three parameters of the models 
in the pressureless dust and electromagnetic radiation 
cases (w = 0, 1/3). 



In this connection, we would like to point out that, 
concerning the cases A ~ (a/a) 2 and A ~ p, it was 
already mentioned by Vishwakarma [16] that the esti- 
mates of the parameters for flat models are the same. 
Therefore, in view of this, the main purpose of the 
present paper is to reexamine the status of the phe- 
nomenological approach of the dynamical A -term and 
to provide more general result by including one more 
case A ~ a/ a into a systematic analysis. However, 
though there are innumerable phenomcnological A- 
decay laws available in the literature (see for exhaustive 
lists [3,17]), this particular case, viz. A ~ a/a, is not 
included there. This case has been so far, separately, 
taken up by Arbab [13-15] and also by Overduin and 
Cooperstock [3] with a different approach. We have, 
among other candidates of the list, purposely omitted 
the popular cases like A <~ tr 2 and A <~ a~ 2 since 
the first one exactly coincides with that of the case 
A ~ (a/a) 2 as t <~ H^ 1 where H is the Hubble param- 
eter, which is defined as H = a/ a and was extensively 
studied by several authors [14-16,18-28]. The second 
case, A ~ a~ 2 , which was first suggested through di- 
mensional arguments related to quantum cosmology by 
Chen and Wu [29] and also results from a contracted 
Ricci collineation along the fluid-flow vector [30,31], is 
dropped here since this case does not suit for our present 
scheme as will be clear from the field equations of the 
next section. 



Based on all the available observational information, 
some physical features have been explored through the 
cosmological parameters, which are in good agreement 
with the observationally obtained present data of the 
Universe. These results are discussed in Sec. 6 consider- 
ing both the present accelerating and the past deceler- 
ating Universe. Before this, we show the ranges of the 
parameters a, (3 and 7 involved in different A -models 
in Sec. 5, whereas the equivalence of the A-models is es- 
tablished in Sec. 4. Sections 2 and 3 are related to the 
Einstein field equations and their general solutions for 
different A-dependent models. In the concluding Sec. 7, 
some discussion is presented. 



2. Einstein's field equations 

Let us consider the Fricdmann-Lemaitre-Robertson- 
Walker (FLRW) metric 



ds 2 = -dt 2 + a{tf 



dr 2 



kr 2 



+ r 2 (d9 2 + sin 2 6dcj> 2 ) 



(1) 



where the curvature constant k = — 1,0, +1 for hyper- 
bolic, flat and closed models of the Universe, respec- 
tively. 

The Einstein field equations are given by 



2 y 



-8irG 



A 



(2) 



where A is the so-called cosmological constant, assumed 
here to be time-dependent, viz., A = A(£), and c, the 
velocity of light in vacuum, is assumed to be unity (we 
thus use relativistic units). 

For the spherically symmetric metric considered 
above, the Einstein field equations with a time-dependent 
cosmological constant yield the following two equations, 
called the Friedmann equation and the Raychaudhuri 
equation: 
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The energy conservation law can be written as 



8irG{p + p)- 



8ttG 



A 
3"' 



(3) 
(4) 

(5) 



Let us choose the barotropic equation of state 

p = wp (6) 

where the parameter w can take the constant values , 
1/3, —1 and +1 for dust, radiation, vacuum fluid and 
stiff fluid, respectively. 

Using Eq. (6), Eq. (4) transforms to 



a 4ttG.., o , A 
a + — (! + *«>/>= 3" 



(7) 



Differentiating Eq. (3) with respect to the time coordi- 
nate t and using Eqs. (4)-(7) to eliminate p, we finally 
obtain the following equation: 



1 



w 



1 + 3w 



1 



k 



w 



1 + 3w 



A. 



(8) 



This is the dynamical equation relating the cosmic scale 
factor a to a known value of the dynamic cosmological 
term A . It can readily be observed from the above equa- 
tions (7) and (8) that A depends on the factors a/a, p, 
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(a/a) 2 and a 2 in a specific way. However, the infla- 
tion theory of the Universe predicts and the CMB detec- 
tors such as BOOMERANG [32-34], MAXIMA [35-37], 
DASI [38], CBI [39] and WMAP [40,41] confirm that the 
Universe is spatially flat. Therefore, the A ~ a~ 2 case, 
which is not suitable for k — 0, is omitted here. How- 
ever, for a detail study of this case, viz. A ~ a~ 2 , inter- 
ested persons may consult the works done by Chen and 
Wu [29], Abdussattar and Vishwakarma [30] and Vish- 
wakarma [31], mentioned earlier, as well as [4,5,42-48]. 
We shall therefore consider the phenomenological mod- 
els related to the cases (a/a) 2 , d/a and p only for A 
and try to find solutions which will help us to model and 
explore the features of the Universe. 



3. Cosmological models for an 
accelerating Universe 

If we use A = 3a(d/a) 2 = 3aH 2 , where a is a con- 
stant and H is the Hubble parameter, then for the flat 
universe ( k = ) Eq. (8) reduces to 



2ad + (1 + 3w - 3wa - 3a)d 2 = 0. 

Solving Eq. (9), we get our general solution as 

a(t) = Cl t 2 /^- a ^ +w \ 
Pit) 



A(t) 



6ttG(1 - a)(l + w) 2 
4a _. 
3(l-a) 2 (l + w) 2 * ' 



(9) 

(10) 
(11) 
(12) 



where C\ is an integration constant. 

It is evident from Eqs. (10)— (12) that a ^ 1 for phys- 
ical validity. Moreover, a repulsive A demands positive 
a via Eq. (12) while Eq. (11) shows that, for positive 
p, the parameter a must be less than 1 and imposes 
the constraint < a < 1 . The case a > 1 is either 
nonphysical or incompatible with a time-dependent A. 
This is because a solution with a variable A is possible in 
the presence of matter only when T n - ^ [49]. Again, 
since we are dealing with a non-zero A , we have a ^ . 
This means that, when A <~ (d/a) 2 , the expansion of 
the Universe never stops as long as A ^ . 

Similarly, if we set A = (3(d/a) and A = 8nGjp, 
where (3 and 7 are free parameters, then, for k = 0, 
it can be very easily shown that in both cases the scale 
factor follows the same type of power laws as in Eq. (10) 
while, just as in Eqs. (11) and (12), p(t) and A(t) are 
inversely proportional to t. It may be mentioned that 
for physical validity either (3 < or (3 > 3 for the 
present model. On the other hand, for a non-negative, 
repulsive A one needs to impose the condition 7 > 
while for positive p it should be 7 > — 1 . This means 
that 7 is always a positive quantity. 



4. Equivalence of three forms of 
dynamic A 

Now, let us explore the interrelations between a , (3 and 
7 and hence the equivalence of different forms of the 
dynamic cosmological terms, viz., A <~ (d/a) 2 , A ~ d/a 
and A ~ p. 

From Eq. (10), differentiating it and then dividing by 
a, we get 



t = 



3(1 - a)(l +w)H 



(13) 



where H is the Hubble parameter, as mentioned ear- 
lier, and hence, for specific values of a and w, Eq. (13) 
shows that H <~ i -1 . This point was indicated in the 
Introduction, and therefore we have omitted the case 
A <~ t~ 2 from the present investigation. 

Using Eq. (13) in (11) and the definition of the cos- 
mic matter density parameter ^ m (— 8ttGp/3H 2 ) , one 
gets 



1 



a 



(14) 



where Q ma is the cosmic energy density parameter for 
the a-related dynamic A-model. 

Again, using Eq. (13) in (12) and the definition of the 
cosmic vacuum energy density parameter = A/3H 2 , 
we have 



n 



Aa 



a, 



(15) 



where, in a similar fashion, il\ a is the vacuum energy 
density parameter for the a-related dynamic A-model. 
Addition of Eqs. (14) and (15) yields 



^ma ~t~ ^Aa 1 



(16) 



which is the relation between the cosmic matter- and 
vacuum density parameters for a flat (k = 0) universe. 

Equations similar to (16) can be obtained for (3- and 
7-related models as well. Thus, without loss of general- 
ity, we can write 



^Aq — ^A/3 — — 



(17) 
(18) 



where Q m and are the cosmic matter and vacuum 
density parameters. Therefore, in the absence of curva- 
ture, one can obtain the general relation 



n = Q m + f2y\ = i. 



(19) 



This analytical result is consistent with the observa- 
tional constraint on the total energy density f2 of the 
Universe, where ft = LOOto' 2 ^ due to the MAXIMA-I 
flight and COBE-DMR experiment [37], fl = 1.05±0.08 
obtained from CBI-DMR observations [39], and il — 
1.01 ± 0.03 (68% CL) measured from the first acoustic 
peak in the angular power spectrum of CMB fluctua- 
tions [50]. 
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Now, Eqs. (17) and (18) enable us to interrelate a, 
(3 and 7 with Q m and £l\ as 



a 

a = 
7 = 



6n A 



20 A - O m (l + 3w)' 



fir. 



(20) 
(21) 

(22) 



This result for the dust case of Eq. (21) corresponds to 
Arbab's [14]. Thus, we find that the free parameter a 
here is nothing but the cosmic vacuum density param- 
eter whereas 7 is the ratio of the cosmic vacuum and 
matter density parameters which, by virtue of Eq. (19), 
provides 



n = (i + 7 )o ro = i, 



(23) 



which is another relation for the total cosmic energy 
density in the case of a flat universe. 

All the above general relations for a, (3 and 7 in 
terms of f2 m and tt\ also hold for the particular cases 
of dust (w — 0) and radiation (w = 1/3). It is inter- 
esting to note that, while the relations of a and 7 with 
the cosmic matter and vacuum density parameters are 
independent of w , the relations of (3 with Q m and Oa 
are w -dependent. 

It can easily be shown that the particular solutions of 
A <~ (a/a) 2 model for dust and radiation cases become 
identical with their corresponding counterparts for the 
other models in terms of the time dependences of a , p 
and A when expressed in terms of fl m and fl A . There- 
fore, these results imply that in £l m and fl A there are 
no distinctive features which could distinguish between 
the different forms of dynamic cosmological models, viz., 
A ~ (a/a) 2 , A ~ a/a and A ~ p. Thus, starting from 
any of our A-models, since they are equivalent, we can 
arrive at the other relations. 

Now, from Eqs. (20)-(22), we find that the parame- 
ters involved in the three dynamical relations are con- 
nected by 



a 



7 



(3(1 + 3w) _ 
3(/3i(; + 0-2) ~ I+7' 



(24) 



This again shows that the three forms A = 3a(d/a) 2 , 
A = (3(a/a) and A = 8irG-fp are equivalent, and the 
three parameters a, (3 and 7 are connected by the re- 
lation (24). Thus it is possible to find out the identical 
physical features of others if any of the phenomenologi- 
cal A relations is known. It can easily be seen that, for 
the dust case (w — 0), Eq. (24) relates a and (3 as 



a = 



308-2)' 



(25) 



which is Arbab's result [15]. Moreover, it can be ob- 
served that our 7 is identical to Majernik's k [51,52], 
where 

! 1 = 2^ (26) 



for the present situation in view of Eq. (22). He has also 
shown that this result, Eq. (26), is derivable from an 
ansatz by which A is proportional to the stress-energy 
scalar T = , the trace of the stress-energy tensor of 
ordinary matter T- , and is Lorentz- invariant. In this 
regard, following Majernik [51], it can be mentioned 
here that determination of the parameter 7 entirely de- 
pends on the cosmic matter density parameter or both 
the matter and vacuum density parameters. Thus this 
relation constrains the value of 7 and will be discussed 
in the next section. 

Another point to be mentioned here that Eq. (19) 
and hence (23), via (22), is nothing but another form of 
the Friedmann equation (3) for the flat Universe. Thus, 
it is interesting to note that Eq. (26) also represents the 
Friedmann equation. Therefore, starting from any of 
our A-models, since they are equivalent, we can arrive 
at the Friedmann field equation without any special as- 
sumption. 



5. Ranges of the parameters a, (3, 7 

Recent measurements have given a wide range of val- 
ues for Q m0 and Q A0 . The first flight of the MAXIMA 
balloon-borne experiment (MAXIMA-I) combined with 
COBE-DMR resulted in 0.25 < O m0 < 0.50 and 0.45 < 
Oao < 0.75 [37]. Observations of SNela combined with 
the total energy density constraints from CMB [50] and 
combined gravitational lens and stellar dynamical anal- 
ysis [53] lead to f2 m o ~ 0.3 and Oao ~ 0.7. The pin- 
point values of these parameters as obtained by Siev- 
ers et al. [39] and Spergcl ct al. [54] are [O m o, £!ao] 
[0.34 ±0.12,0.67^0.13] and [O m0 , O\o]= [0.249 ~ 



+0.024 
0.031 > 



0.7191 

0.029 L respectively. These and other results are 
listed in Table 1. 

Considering the values in Table 1, we particularly 
prefer the matter density parameter as £l m o = 0.330 ± 
0.035 [49,50,55,56]. This gives us an opportunity to 
obtain ranges of ao , A) an d 70 (the values of a , (3 and 
7 at the present epoch) which can, using Eq. (24), be 
obtained as 0.635 < a a < 0.705, 3.417 < (3 a < 4.674 
and 1.739 < 70 < 2.389 in the dust case. Thus we find 
that using our models we are able to obtain the range 
of (3 smaller than Arbab's [15] which was 3 < (3 < 4.5 
for dust. 

Again, if we recall the quintessence equation of state 
PQ = W QPQ where wq — — Oa , we can easily obtain the 
relations between (3 and wq as (3 — 6wq/(1 + 3wq) for 
w = . Using the above range of (3$ , we can calculate 
the range of wq in dust as —0.705 < wq < —0.635. It 
is interesting to note that the above result is consistent 
with the accepted range of wq which is — 1 < wq < . 
However, in the present investigation we arc not con- 
cerned with the quintessence case and show the range 
of wq as a check only. 
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Table 1: Values of il m0 and f^Ao from various observa- 
tional sources 



Source & Reference 


Year 


!!m0 




rLtstatmou et al. [oi\ 


tyya 


fl nr + 0.18 
U - Zd -0.12 


0.63I 023 


(blNela + UVLBJ 








Riess et al. [1] 


1998 


n 0,1+0.56 

0.24I 024 


72+ ' 72 










Pcrlmutter et al. [58] 


1999 


0.4 ± 0.1 


0.7 


^blNela + Ijob j 








Balbi et al. [37J 


2001 


0.25—0.50 


0.45 — 0.75 


(\/\ A VTT\/T A T I Pn'Dl? 1 ! 








IXCfJUltJ [OU\ 




u.ou 




(SNela + CMB) 








Koopmans et al. [53] 


2003 


0.30 


0.70 


(Lens + SD) 








Sievers et al. [39] 


2003 


0.34 ±0.12 


0.6712:1° 


(CBI) 






Barris et al. [59] 


2004 


0.33 


0.67 


(IfA Deep Survey) 








Astier et al. [60] 


2006 


0.31 ±0.21 


0.80 ±0.31 


(SNLS) 








Spergel et al. [54] 


2006 


0.2491°;^ 


o.7i9lg:SJ 



(WMAP + SNLS) 



6. Features of the models 

6.1. Physical parameters at the present 
accelerating epoch 

Now, a search for the status of A rests on some obser- 
vational results from high-rcdshift type la Supernovae 
(SNela), the cosmic microwave background radiation 
(CMBR) and other sources which inform us that the 
present Universe is composed of about 30% of ordinary 
matter and 70% of dark energy. Thus, in the present 
Universe, the vacuum density parameter Oao is domi- 
nant over the matter density parameter il m0 . Determi- 
nation of the Hubble parameter, a measure of the rate 
of cosmic expansion, has been done by several authors 
based on different values of density parameters as shown 
in Table 1. However, it is to be noted that there exists 
a certain amount of uncertainty in the value of Ho , as 
is obvious from Table 2. 

The data of Table 2 indicate that the present value of 
the Hubble parameter is, in general, centralized at 72±8 
kms _1 Mpc _1 . Even the most recent value (73.4l 3 'g) 
as obtained from WMAP by Spergel et al. [54] lies 
well within this range. Assuming this value of Ho and 
O ro0 = 0.330 ± 0.035 [49,50,55,56], the present age 
(to), the present matter density (po), the present value 
of the cosmological term (Ao) and the value of the de- 
celeration parameter at the present era ( qo ) have been 
calculated using our equivalent models. All values of p 
and go arc m nice agreement with the modern concept 
of an open, accelerating universe. Moreover, the values 
of A support the idea of a small non-zero cosmologi- 
cal parameter which is slowly decreasing in time, and 



Table 2: Values of H and to from various observational 
sources 



Source & Reference 


Year 


km/(s- Mpc) 


to 

Gyr 


Sandage, Tammann [63] 


1997 


55 ± 5 


13.5 ± 2.5 


Sandage et al. [64] 


1998 




I3.5I3 


Parodi et al. [65] 


2000 


58 ±6 




Birkinshaw [66] 


1999 


60 ± 10 




Salia et al. [67] 


1999 


60 ±2 




Jha et al. [68] 


1999 


64+* 

w — 6 


14.1 ± 1.6 


Perlmutter et al. [2] 


1998 


65 


14.9 


Riess et al [1] 


1998 


65.2 ± 1.3 


14.2 ± 1.7 


Tonry et al. [691 

J L J 


2003 


58.8 - 72.3 


13.1 


Sievers et al. [39] 


2003 


66 ± 11 


14.2 ± 1.3 


Tegmark et al. [70] 


2003 


66+«i 

—6.4 


14.1+12 

— 0.9 


Rebolo [50] 


2003 


60 - 75 


13.0 ± 1.5 


Altavilla et al. [71]] 


2004 


(68-74) ± 7 




Cardone et al. [72] 


2004 


56-88 


13.1-14.3 


Freedman et al. [73] 


2001 


72 ±8 


13 


Knox et al. [74] 


2001 




14.0 ± 0.5 


Ferreras et al. [75] 


2001 




13.2+H 

^ — 0.8 


Alcaniz [56] 


2004 




13.7 ±0.2 


Freedman, Turner [76] 


2004 


72 ±7 


13.0 ± 1.5 


Spergel et al. [41] 


2003 


72 ±5 


13.4 ± 0.3 


Spergel et al. [54] 


2006 


73.4+oi 

1 —3.8 


13.73 + nl? 
^ ~" — . 1 7 


Freedman [77] 


1996 


73 ± 6 




Riess et al. [78] 


2005 


73 ±9 




Peebles, Ratra [79] 


2003 


73 ± 10 




Blakeslee et al. [80] 


1999 


74 ±4 




Koopmans et al. [53] 


2003 


751^ 




Efstathiou [81] 


1995 


80 




Jacoby et al. [82] 


1992 


80 ± 11 




Freedman et al. [83] 


1994 


80 ± 17 




Pierce et al. [84] 


1994 


87 ±7 





at present A lies within 1 x 10~ 35 s~ 2 — 2xl0~ 35 
s~ 2 , which agrees with the results of Carmeli [61] and 
Carmeli and Kuzmenko [62], where they obtain the 
value of 1.934 x 10~ 35 s~ 2 . All values of p are one 
order of magnitude smaller than 10~ 29 g/cm 3 , the crit- 
ical density of the Universe. For the matter-dominated 
case, various results can be obtained for t by finding 
H for different Q m0 (see Table 3 for detail). 

It is clear from Table 3 that, for the lower value 
of f2 m0 , the age becomes very high, whereas a higher 
value of the matter density parameter, say, f2 m o = 0.46 
provides a more realistic result for the age of the Uni- 
verse with gradual increase of the Hubble parameter. 
The best result is therefore obtained for Q m0 = 0.46 
(the upper limit of Sievers et al. [39]) and H = 80 
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Table 3: Age of the Universe from the present models 



( ) 


rj 


to 




krns/ (s-Mpc) 




0.33 


64 


27.79 




72 


24.70 




ou 


99 9^ 


0.365 


64 


25.13 




72 


22.33 




80 


20.10 


0.40 


64 


22.93 




72 


20.38 




80 


18.34 


0.46 


64 


19.94 




72 


17.72 




80 


15.95 



kms^Mpc" 1 (the upper limit of Refs. [56,69,74,75]) 
is 15.95 Gyr. This result exactly coincides with the up- 
per limit of the value of Ricss et al. [1] which is 14.2±1.7 
Gyr as obtained from SNela observations and also very 
close to the values obtained by Sievers et al. [39] and 
Tcgmark et al. [70] as predicted by the WMAP data 
and CMB observations. 

6.2. Physical parameters in the past 
decelerating period 

Using Eq. (10), one can obtain an expression for the de- 
celeration parameter q as 



aa 



3(1 - a)(l+w) 



- 1. 



Thus, for an accelerating universe 
1 + 3w 

a > 



3(1 + w) 



(27) 



(28) 



From Eq. (28), it is evident that for dust (w = 0) an 
accelerating universe requires a > 1/3. Now, a be- 
ing the cosmic vacuum density parameter by virtue of 
Eq. (20), we find that our model fits an accelerating uni- 
verse since the modern accepted value of 12ao is about 
0.7 (Refs. [72,85-87] and see also Tabic 1) and is much 
larger than 1/3. Thus Eq. (19) implies that the value 
of fl m is 0.3, which provides qo = —0.50 ± 0.05 for 
the dust case which can nicely accommodate the cur- 
rently accepted value related to the accelerating Uni- 
verse [57,88,89]. Again, q will be positive if a is less 
than 0.3. Thus, for a decelerating universe, the cosmic 
vacuum density parameter should be smaller than 0.3, 
which is also consistent with the modern ideas. There- 
fore, we find that, within our models, one can investigate 
accelerating as well as decelerating phases of the cosmic 
expansion since q depends on a . 

Now, it has already been mentioned that the expand- 
ing Universe, which is about 14 Gyr old, entered into 



the present accelerating phase about 7 Gyr ago. There- 
fore, for about 8 Gyr earlier, i.e., when the Universe 
was about 6 Gyr old, it was passing through a period 
of deceleration. Let us try to estimate the values of 
some of the physical parameters when the age of the 
Universe was 6 Gyr. From Eq. (27) it is easy to obtain 
q = (1.5f2 m — 1) for w = 0. We have already seen that 
q will be positive for TO > 0.66. Putting TO = 0.67, 
we find that q > 0, i.e. the Universe was indeed in a 
decelerating phase. Assuming Q, m = 0.67 and t = 6 
Gyr, we can estimate the values of H , p and A, which 
are given by H ~ 179 km s _1 Mpc _1 , p = 3.3 x 10~ 29 
g cm" 3 and A = 2.74 x 10~ 35 s~ 2 , respectively. Sim- 
ilarly, for the radiation case, f2 m and q are related by 
q = (2Q m — 1). Thus q will be positive for fl m > 0.5. 
If we assume Q TO = 0.6 and t = 6 Gyr, then the esti- 
mated values of H, p and A are 135 km s _1 Mpc _1 , 
2.08 x 10~ 29 g cm- 3 and 2.33 x 10" 35 s" 2 , respec- 
tively. We thus find that in both cases (to = 0, 1/3) p 
was above the critical density, which means that the ex- 
panding Universe with a decelerating mode had a closed 
geometry. The value of the Hubble parameter, a mea- 
sure of the expansion rate of the Universe was slower 
in the radiation era than that of matter-dominated era. 
Also, we find that the value of A was slightly above its 
present value, which justifies the idea of a dynamic A 
which decreases very slowly with time. Finally, assum- 
ing the present value of the Hubble parameter as 72 km 
s Mpc , we see that the rate of decrease of H is 
about 13 km s _1 Mpc -1 Gyr -1 for w = and about 8 
km s _ 1 Mpc ~ 1 Gyr ~~ 1 for w = 1/3. 

7. Discussion 

In the present investigation, choosing some specific 
forms of dynamical A , we were able to show the equiv- 
alence of those forms in terms of the solutions obtained. 
While Arbab [14] has shown the equivalence of the same 
three A -models in the context of a built-in cosmologi- 
cal constant of Rastall [18] and Al-Rawaf-Taha [26] type 
models of modified general relativity, we have shown the 
equivalence of the models with respect to their charac- 
teristic solutions in the framework of Einstein's general 
relativity. It has already been mentioned that A ~ H 2 
and A ~ t~ 2 are identical. In this context, it is inter- 
esting to note that since d/a is equal to H + H 2 which 
is again ~ t~ 2 , then A ~ a/a is also identical with the 
above-mentioned cases. This is reflected in our solution 
sets via Eq. (24). Moreover, since a/ a = H + H 2 , the 
A ~ d/a model can be thought of as a combination 
of two models, viz., A ~ H and A ~ H 2 . Thus the 
A ~ d/a and A ~ a 2 /a models become identical when 
H = 0. Now, H = implies a constant H, which in 
turn implies exponential expansion and hence an infla- 
tionary scenario. Thus the idea of inflation is inherent 
in the phenomenological model A ~ d/a. Moreover, the 
A ~ a 2 /a and A ~ d/a models cannot exist as separate 
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entities during inflation. 

We have also established a relation between a, j3 
and 7, the three parameters of the three forms of A 
which ultimately yields Q m + Qa = 1 , the relation be- 
tween the cosmic matter and vacuum density param- 
eters for a flat universe. It can be shown that this 
particular relation between the density parameters also 
holds in the radiation case. On the other hand, since 

= Pa/Pc and Q m = p m j p c , it is clear from Eqs. (20), 
(21) and (22) that (3 and 7 are independent of p c , the 
critical density of the Universe, whereas a depends on 
the critical density. It can also be shown that, while (3 
and 7 decrease with the age of the Universe, a increases 
as time passes. 

Moreover, the present models represent a flat, ac- 
celerating Universe and do not suffer from the low- 
age problem like many FRW models. Also, since the 
present Universe is dark energy dominated, and the 
closest approximation to t ~ 20.10 Gyr for the matter- 
dominated case can be obtained for the specific choice of 
O ro0 = 0.330 + 0.035 and H = 72 + 8 km s^Mpc" 1 , 
our models point at the upper accepted limits of fl m o 
and Ho ■ As has been shown earlier in Table 2 that 
there is a certain amount of uncertainty in the value 
of H , the lower bound being 50 km s _1 Mpc _1 (see 
also Table 2 of Ref. [64] for some more cases of a lower 
bound), whereas the upper bound is 97 km s _1 Mpc -1 . 
The related values of the age of the Universe for our 
models with these two extreme Hubble parameters are 
25.52 Gyr and 13.15 Gyr, respectively, when the mat- 
ter density parameter is 0.46. Therefore, our calculated 
value of to , which seems to be a bit over-aged and also 
favours a A -dominated universe, can be accommodated 
to the accepted age of the Universe within the error 
bar. Even though the values of t in case w = show 
an excess with respect to the presently accepted age of 
the Universe but, anyway, there is no low-age problem. 
Indeed, these values are much higher than the age of 
the globular clusters which is 12.5 ± 1.2 Gyr [90-95]. 
In this connection, we can also note that examples of 
higher age are not unavailable in the contemporary lit- 
erature [49,96-102]. For example, Vishwakarma [49] 
obtained, for O ro0 = 0.330 ± 0.035 and H = 72 ± 7 
km s _1 Mpc _1 , a remarkably high age of the Universe, 
t w 27.4 ± 5.6 Gyr! But it is evident from Table 2 
that, whatever be the values of the Hubble parameter, 
the experimental results for the age of the Universe lie 
around 14 Gyr. For the present phenomenological mod- 
els (including Vishwakarma's case [49]), the age of the 
Universe is inversely proportional to the Hubble param- 
eter. This provides a reasonable age of the Universe 
only for a higher value of the Hubble parameter, which 
is also clear from the above discussion. This is obviously 
a drawback of the present models unless a higher value 
of Hq is observationally established in the future. 

In this regard, we would like to discuss the causal 



connection of our models. We know that the proper 
distance L(t) to the horizon is given by 

L(t) = a(t) [ dr/a(r), 

and if the integral diverges, the model is causally con- 
nected. Since, for A ~ (d/a) 2 , the scale factor a(t) is 
given by Eq. (10), the proper distance L(t) diverges if 
1/3 < a < 1 or, in other words, if 1/3 < Qa < 1 for 
the dust case and a < 0.5, i.e., fl\ < 0.5 for radiation. 
Since the present Universe is matter-dominated and the 
observational results indicate that at present J7a ~ 0.7, 
the Universe is causally connected in our A ~ (d/a) 2 
model. Besides, since it has already been shown that 
the present three phenomenological A-models are equiv- 
alent, the causal connection of the Universe indicated in 
the above model implies that the other models are also 
causally connected. 

It should be mentioned that Arbab [14] has put 
his models to the neoclassical tests like luminosity dis- 
tance, angular diameter distance and gravitational lcns- 
ing, whereas we have tested the viability of our models 
through age determination and some other measure- 
ments. So, in that respect our work can be thought of 
as complementary to Arbab's investigation [14]. Per- 
spectively, we are studying new forms of a dynamic 
cosmological term, such as those from the renormal- 
ization group, and its confrontation with astrophysical 
observational data sets [103,104], expectating that this 
global description can help one to better understand 
the mysterious dark energy nature and to alleviate the 
long-standing cosmological constant problem. 

Finally, it is to be noted here that, in general, to is a 
function of time [79,105,106]. But the current observa- 
tional data can hardly distinguish between time- varying 
and constant equations of state [107, 108], as demon- 
strated in some works ( [108] and references therein). 
For this reason, w is usually assigned a constant value, 
as has been done by Caldwell et al. [109] while deal- 
ing with a relation between scalar field models and the 
XCDM parameterization. Likewise, in the present work 
related to phenomenological A models, without show- 
ing a complete time evolution of w (which is no doubt 
a better representation), some specific cases are high- 
lighted corresponding to constant w . A more accurate 
analysis may be made in a later work by considering 
w = w(t) . 

Appendix 

A comparative analysis of models 

A = a{a/af, A = f3(a/a), A = jp. (A.l) 

can be, alternatively, done as follows. It can be shown 
that all three models do not differ from one another from 
both mathematical and physical points of view. 
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Let us denote 
h = a{a/a) 2 , 



fa = P(a/a), 



fa = IP- (A.2) 



Then Einstein's equations for spatially flat space-time 
used in this work (Eqs. (3) and (7)) can be written in 
the form 



h + Ofa-^fa 



o/ 1 + / 2 + ^d 



3 



0. 



I A = o. 



(A.3) 
(A.4) 



This is a set of linear algebraic equations for fa , fa and 
fa with constant coefficients. If we add to this set any 
linear equation of the form 



a fa + bfa + cfa = dA 



(A.5) 



with constant coefficients a, 6, c, d, then the Einstein 
equations and this new equation make a closed set of 
three linear algebraic equations for fa, fa, fa . Its 
solution for any real a , b , c , d can be written as 



h 



i 



-A, 



1 



fa 



1 

1 



-A. 



(A.6) 



A comparison of Eqs. (A.2) and (A.6) clearly shows 
that they are equivalent to Eq. (A.l). However, an im- 
portant point is the arbitrariness of a, b, c and d. 
Therefore, using one of the relations in Eq. (A.l) au- 
tomatically leads to one of the other two relations, and 
consequently there is no difference in the dynamic be- 
haviour of these models. 
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